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We use the recently derived density of states for a particle confined to a spherical well in three 
dimensional fuzzy space to compute the thermodynamics of a gas of non-interacting fermions con¬ 
fined to such a well. Special emphasis is placed on non-commutative effects and in particular 
non-commutative corrections to the thermodynamics at low densities and temperatures are com¬ 
puted where the non-relativistic approximation used here is valid. Non-commutative effects at high 
densities are also identified, the most prominent being the existence of a minimal volume at which 
the gas becomes incompressible. The latter is closely related to a low/high density duality exhib¬ 
ited by these systems, which in turn is a manifestation of an infra-red/ultra violet duality in the 
single particle spectrum. Both non-rotating and slowly rotating gasses are studied. Approximations 
are benchmarked against exact numerical computations for the non-rotating case and several other 
properties of the gas are demonstrated with numerical computations. Finally, a non-commutative 
gas confined by gravity is studied and several novel features regarding the mass-radius relation, 
density and entropy are highlighted. 

PACS numbers: ll.lO.Nx 


I. INTRODUCTION 

The structure of space-time at short length scales is probably the most challenging problem facing modern physics 
[T]. Although there seems to be growing consensus that our notion of space-time at the Planck length needs major 
revision, there seems to be much less consensus on the appropriate description of space-time at these length scales. 
One scenario, originally proposed by Snyder [2], is that of non-commutative space-time. Although his motivation 
for this proposal was an attempt to avoid the infinities encountered in quantum held theories, the notion of non- 
commutative space-time has recently been revived by the rather compelling arguments of Doplicher et al [3] and the 
emergence of non-commutative coordinates in the low energy limit of certain string theories [3]. 

Since the energies required to probe the Planck length scale only occurred in the very early universe, we do not 
have any direct empirical data to guide our search for the structure of space-time at the Planck length scale. Possible 
sources of indirect data are the cosmic micro-wave background radiation (GMB), other cosmological phenomena such 
as dark matter and dark energy, high energy cosmic radiation and cold, dense astro-physical objects such as neutron 
stars and white dwarfs. In particular it has been observed before [S] that the thermodynamic behaviour of two 
dimensional Fermi gases is drastically altered at high densities and low temperatures. As these computations were 
performed for two dimensional systems, not much can be inferred about the behaviour of such real physical objects. 

This motivates the current paper in which the thermodynamic behaviour of Fermi gases in three dimensional 
non-commutative space is investigated. To avoid the difficulties related to the breaking of rotational symmetry 
when the simplest constant commutation relations are assumed for the coordinates [6], we rather adopt here fuzzy 
sphere commutation relations for the coordinates. The recent solution for the spectrum of the non-commutative 
three dimensional fuzzy well [7] is then used to compute the thermodynamics of these gases, confined to a finite 
volume, in three dimensions. Particular attention is paid to the behaviour of these gases under conditions where 
non-commutative effects are expected to be relevant, i.e. at high temperatures and in dense, cold Fermi gases. These 
computations, based on solving the non-commutative Schrodinger equation, are still non-relativistic and certainly 
not realistic for the ultra-relativistic conditions prevailing at high temperatures and in dense astro-physical objects. 
However, although the non-relativistic description may affect the quantitative behaviour, one does not expect it to 
alter the qualitative behaviour, such as incompressibility. We therefore expect to learn a considerable amount about 
the qualitative behaviour of such systems in a non-relativistic description, which is considerably simpler than the 
corresponding more realistic, but also much more difficult, relativistic description, which is as yet still lacking. We do, 
however, also compute the non-commutative corrections at low densities and temperatures where the non-relativistic 
approximation is certainly appropriate. These corrections may offer useful guidance to attempts to observe these 
phenomena at low energies and densities. 

In this paper the focus falls on Fermi gases, where one expects non-commutative effects to be most prominent. 
One can, of course, also study the thermodynamics of Bose gases, but as one expects non-commutative phenomena 
to be less prominent there, we postpone this study to a future paper. Furthermore, at low densities these gases are 
essentially classical or Boltzmann, where quantum and statistical effects are unimportant, and the non-commutative 
corrections computed in this limit should show up regardless of whether we consider a Bose or Fermi gas. 
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The paper is organised as follows: Section [T^ reviews the basic formulation of quantum mechanics in fuzzy space 
and the solutions of the infinite fuzzy well. Section III sets up the formalism to study the thermodynamics of the 
non-commutative Fermi gas. In particular it is demonstrated that a remarkable infra-red/ultraviolet duality in the 
spectrum results in a similar duality between the low and high density limits of a non-commutative Fermi gas. Section 


IV focuses on the computation of the central thermodynamic quantities, firstly, in the low density, low temperature 


as well as low density, high temperature limits. The low/high density duality is then applied to infer from this the 
behaviour of these quantities in the high density limit. Both non-rotating and slowly rotating gases are discussed. 
Section [V] discusses a Fermi gas confined by gravity and in particular the effect of non-commutativity on the mass- 
radius relation, density and entropy. Section [Vl| summarises our results and draws conclusions. 


II. SPECTRUM OF THE INFINITE FUZZY WELL 
A. Basic formalism 

We begin by summarizing the formalism of quantum mechanics in non-commutative three-dimensional space, as 
set out in mu- The starting point is to impose on the spatial coordinates the su{2) commutation relations 

[X,,Xj] = 2iee,,kXk ( 1 ) 

where 9 is the non-commutative length parameter. The Casimir operator X'^ = XiXi is naturally associated with the 
radial distance squared and its eigenvalues are determined by the su(2) representation under consideration. 

The first step in setting up the non-commutative quantum system is to build a representation of the non- 
commutative coordinate algebra 0 on a Hilbert space "He, referred to as the configuration space. In this case 
we want to mimic as a collection of spheres (an onion structure) where the squared radii of the spheres are quan¬ 
tised according to the eigenvalues of the su{2) Casimir. For this He must carry a single copy of each su(2) irrep. A 
convenient, concrete realisation of this configuration space is provided by the Schwinger construction, which utilises 
two sets of boson creation and annihilation operators to build a representation of su{2): 

[aa,a^p]= Sap and [cq,, a^] = [a)),, aj,] = 0 with a,/3 = I,2. (2) 

The corresponding boson Fock space is identified with He- The coordinate operators are realised as 

A, = Oa\a^^lap (3) 

where } are the Pauli spin matrices. The Casimir operator reads X"^ = XiXi = 9'^N{N + 2) with N = a|ai - 1-0202 
from which it is clear that each su(2) representation occurs precisely once in He- As a measure of radial distance we 
will use 


f = 9{N + I). 


( 4 ) 


The quantum Hilbert space Hq is now defined as the algebra of operators generated by the coordinates, i.e. the 
operators acting on He that commute with and have a finite norm with respect to a weighted Hilbert-Schmidt 
inner product [7]: 


'^Q = \'^= X! Cj^^-^rn\j,m'){jm\ : tre -h < 00 


( 5 ) 


Here trq denotes the trace over configuration space. 

The quantum mechanical angular momentum operators are realised on Hq as 


^*■0 = with [Li,Lj\ = ihSijkLk. 


The simultaneous eigenstates of and L 3 are of the form 






mi! m2! 


: RiN) : 


nx\ni\ 


( 6 ) 


( 7 ) 
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with j = 0,1, 2,... and m = —j ,..., +j. The summation above is over the non-negative integers (rni,m 2 , ni, 712 ) 
which satisfy mi + m 2 = ni + n 2 = j and mi — m 2 — n-i + n 2 = 2m. The state’s radial dependence is determined by 
the function i?(n) = Cfcn^ through its normal ordered form 


R{N) 


R{N):=Y,Ck 

k=0 


NI 

{N-k)\' 


( 8 ) 


B. Free particle solutions 


In this, and the following two subsections we collect some results relevant to the non-commutative free particle 
and square well problems. Our treatment follows that of section 7 in [7]. We first introduce the non-commutative 
analogue of the Laplacian [8] which acts on a € "Hq as 




Using the form of ipjm in the free particle Schrodinger equation 


2 mi 


-Nelpjrn = 


jm 


(9) 


( 10 ) 


can be reduced to a difference equation for R{n) on the non-negative integers, or equivalently, a differential equation 
for R(n) on n £ [0, 00 ). The two solutions are related by the normal ordering procedure in (|^. Here we only require 
R{n), which is found to be 


where k = 9k, E = and Pn°‘’^\x) is a Jacobi polynomial. This is the non-commutative analogue of the usual 
Bessel function solution of the radial wave equation. By using identity (22.15.1) in [T^ this expression can be shown 
to reduce to the standard commutative result as 0 —>■ 0. An important consequence of 0 is that it places an upper 
bound on the free particle energy E. This can be seen by considering the asymptotic behaviour of R{n) as n —>■ 00 . 
According to theorems (8.21.7) and (8.21.8) in [11] the limit lim„_,.oo R{n) will diverge if the argument of the Jacobi 
polynomial falls outside the interval [—1,1]. For ipjm to be normalisable with respect to the trace inner product in 
1^ K therefore cannot exceed two, which translates into an energy upper bound oi E < ■ 


C. The infinite fuzzy well 

We next consider a spherical well potential in non-commutative space. To this end, consider the operator 

00 n 

Q= 0 '^\k,n-k) {k,n-k\ ( 12 ) 

n—M+1 k—0 

which projects onto the subspace of He spanned by Fock states with total particle number greater or equal to M + 1. 
Through the relation between N and radial distance in Q this amounts to projecting onto configuration space states 
which are localised beyond a physical radius of i? = 6{M -|- 1). For a potential well with height U > 0 and radius R 
the Schrodinger equation now reads 


2mo 


Aotp^ 


•jm 


+ V Q'lpjni — Eljjjjfi. 


(13) 


Through Q this again translates into a difference equation for R{n). In this equation M. = M — j is found to appear 
naturally as the effective radius of the well. The solution for n < A1 J- 1 is just that of the free particle given in ( [TT] ), 
while for n > A1 it is a linear combination of non-commutative versions of the Bessel and Neumann functions. A 
careful analysis of the matching condition dX n = M.,M. + 1 reveals that in the U —?> 00 limit the interior solution R{n) 
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must vanish at n = + 1. According to (111 the bound state energies of the infinite well are therefore determined 

by the zeroes of the Jacobi polynomials. Upon using the symmetry relation Pn°‘’°‘\—x) = (—and the 
freedom to shift the energies by a constant (thereby taking k^/ 2 — 1 to k^/ 2) the quantization condition reads 




(14) 


Note that, due to shifting the energies, the bound states are now arranged symmetrically around E = 0. 


D. Properties of the bound state spectrum 


It follows from (14) that for a specific angular momentum j there are exactly (2j + 1)(M — j + 1) bound states. 
This implies not only a cut-off in energy, but also in angular momentum at jmax = Af. By summing over j the total 
number of single particle states, and therefore the maximum number of spinless fermions the well can accommodate, 
is found to be 




1 A/f3 

= -(M+ 1)(M + 2)(2M + 3) = --b O(M^). 

6 3 


(15) 


As a consistency check on (14), we derive this result in an alternative way. Since the infinite well of radius M 


confines particles within this radius, the configuration space is effectively truncated at states with a radius less 
or equal to M and is thus finite dimensional. The states spanning this finite dimensional subspace of Fock space 
are \ni, 712 ) where ni -b n 2 = = 0,1,... M. Since the quantum Hilbert space is spanned by states of the form 
\N,m'^m) = \j = N/2,m'){j = N/2^m\ where N = 0,1,... M and to,to' = —N/2,... N/2, the quantum Hilbert 
space is also finite dimensional with dimension J2n=o(^ + ^)^’ which agrees precisely with ( |l5| . Thus, counting the 
number of independent solutions of the non-commutative Schrodinger equations matches precisely the dimensionality 
of quantum Hilbert space. 

In order to study the thermodynamics we require the density of states, which is related to the density of zeros of 
the Jacobi polynomial . The latter, in the limit of large j and M, is Md{x,X) where A = j/M and 

[HUS] 


d{x, A) = 


\/l — x'^ — A^ 

TT (1 — x'^ ) 


(16) 


By definition d{x,X) = 0 when \x\ > Vl -b A^. The thermodynamic density of states, still for a fixed j, is then 
p{E,M,j) = {2j + 1)Md{E/E q, j/M)/E q with Eq = and where ^ ranges from —y/\ — A^ to -b\/l — A^. Note 

that the density of states is symmetric under E —>■ This simple fact will have important implications for the 

discussion of the high/low density duality in section (mg Also of interest is the total density of states p(A, M) 
obtained by summing p{E, M,j) over j. This is found to be p{E, M) = NPd{E/E q)/Eq where 


d{x) = / 
Jo 


dX2Xd{x, X) = 


2^1^ 

Stt 


(17) 


Figure illustrates the relevant properties of d{x,X) which are imprinted on the density of states p{E,M,j). The 
cutoffs in the angular momentum ai j = M and energy at E = Eq^/I — {j/M)"^ are clearly visible, as is the symmetry 
of the density of states around zero energy. 


III. THE NON-COMMUTATIVE THREE DIMENSIONAL FERMI GAS 


A. The thermodynamic g-potential 


We now turn to the thermodynamics of a non-interacting gas of spinless fermions which are confined to an infinite 
spherical well with radius R — 9{M + 1) in non-commutative space. We will treat this system in the grand canonical 
ensemble and fix the average value of the z-component of the total angular momentum. The grand canonical potential 
reads 




M 


+i 

E 


log 




j—0 m——j n 


(18) 
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A=jlM 


x—EjEo 


FIG. 1: Figure (a) shows a gradient plot of d{x,X) which defines the density of states p{E, M, j). Figure (b) shows d(x,\) as 
a function of x for the four values of A indicated by dashed lines in figure (a). 


Here {xnj} is the set of M — j + 1 zeros of the Jacobi polynomial (x), as is required by the quantisation 

condition in ( |l4[ ). We have also introduced the non-commutative energy scale Eq = in terms of which we defined 
the dimensionless inverse temperature j5 = Eq j3 and chemical potential p, = The dimensionless angular velocity 
Co' = ^ serves as a Lagrange multiplier for fixing the dimensionless total z-component of the angular momentum 

£ 3 °* = Throughout we follow the convention that quantities denoted by a tilde are dimensionless. 

In the thermodynamic (large M) limit the sums over j and m may be re plac ed by integrals over A = j/M and 
a = m/M while the n summation becomes f dxd{x,X) with d{x,X) given in (16). This produces 


q{M, P, p,Q) = / dX 


r+A 


r^+(A) 


da 


/-A 


e-(A) 


dxd{x, X) log 


1 + e 




(19) 


where Co = MCo' and x±(A) = ±-\/l — A^. This expression serves as the starting point for calculating the various 
thermodynamic quantities. These are obtained from q using 


N=i 

P 


dq 

dp 


^(tot) _ 


T,ui,V 


M 

J 


dq 

dCo 


S 


T,ft,V 


Here S is the entropy, N the number of particles, P = a dimensionless pressure and V = 

dimensionless system volume. The latter is also (to leading order in M) the total number of single particle states 
available to the system and thus represents the maximum number of spinless particles that can be accommodated by 
the system. It is therefore convenient to define 


\dp 


P = 




pv 


( 20 ) 


4 ^ the 




= V = 


M3 

IT 


V 


( 21 ) 


47103 

Note that the physical meaning of A^max is the total number of cells with volume Vq = that fits into the volume 
V. It is therefore also useful and sensible to define a maximum density, which is the density obtained when each cell 
is occupied by exactly one particle, i.e. 


Pmax — T r — 


1 

V)) 


1 

47703 


( 22 ) 


In terms of Pmax the dimensionful density reads 


N 

P ~ ~Y ~ 0niax0 


( 23 ) 
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where p = ^ is the dimensionless density. 

We introduce one further notion that will turn out to be useful later. Suppose we fix the number of particles in the 
system to be N. Then there is a minimum volume that can accommodate this number of particles given by 


N 


14, 


mm 


14iin 

47r03 


(24) 


Here Mmin is the smallest integer larger than or equal to (3iV)^4^ go that the equality between N and t4iin holds to 
0(1/Mmin). In further applications where N and Mmin are assumed to be large, this correction will be ignored. 

It will often be convenient to work in terms of the fugacity z = In terms of this the g-potential of (19) 

reads 


q(M,p,z,w)=M^ 


dX 


+ A 


da 


x+{\) 


dx d{x, A) log 


^-P(x-Lba) 


(25) 


Jo 




Ix-(X) 


where z << I corresponds to the low density limit and z >> I the high density limit. The corresponding thermody¬ 
namic identities are then given by 


N = z 


(dqX 


T (tot) 

^3 


M 

/3 Jf,z,v 


^=q-Nlogz-p(^'] P=X. (26) 


B. The high/low density duality 


Using the symmetry of the density 


of states under x —>■ —x one can trivially rewrite (25) as 


q{M,^,z,6j) 




logz-b g(M, ,0, z \w). 


(27) 


This demonstrates a remarkable duality between the high and low density g-potentials, i.e., we can express the high 
density g-potential (z >> 1) in terms of the low density g-potential (z << 1). This duality, in turn, is a direct 
consequence of the infra-red/ultraviolet duality of the density of states. It is therefore sufficient to compute the 
g-potential in the low density limit z << 1 and infer the high density behaviour from the above duality. 


IV. THERMODYNAMICS OF THE NON-COMMUTATIVE THREE DIMENSIONAL FERMI GAS 


In this section we compute the thermodynamics of the non-commutative Fermi gas. Given the high/low density 
duality above, our strategy will be to first compute the g-potential in the low density limit and then derive the high 
density behaviour. It is also convenient to consider the unconstrained ensemble with w = 0, and therefore = 0, 
and the constrained ensemble with ^ 0 separately. 


A. The unconstrained ensemble with uj = 0 and — 0. 

We start by considering the unconstrained case with a) = 0, which corresponds to an average angular momentum 
= 0. In this case the g-potential reads 


g(M,/3,z) = 


2M3 

Stt 


dx \/\ — x^ log 


1 -b ze 


which follows from combining (16) and ([l7|) in (19). 


(28) 


1. The low density limit. 

In this limit z << 1, but we still have to consider two cases. In the first case the dimensionless chemical potential 
T = j log z < —1 and in the second case 0 > p. = j log z > —1 where, for a low density, (i « —1. 
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We begin by considering A = ^ log z < — 1 in the high and low temperature limits, starting with the high temperature 

limit when /3 << 1. After solving the fugacity to second order in the density, which is required to detect the lowest 
order non-commutative corrections to the pressure, we obtain for the central thermodynamic quantities to the same 
order in the density, lowest order in j3 and after restoring dimensions 


0 = 


P 


1 + 


Pmax 

P = pkT ( 1 + 
S = Nk 


Pmax 

P 

2p 

ma? 


«1, 


1 - log 


p 


2pn 


(29) 


where p is the dimensionful density. Note that the dimensionful temperature here corresponds to kT » Eq so 
that these are ultra high temperatures and not physically realistic. It is therefore also not surprising that non- 
commutative effects feature explicitly as the thermodynamics will be sensitive to the high energy cut-off Eq at these 
high temperatures. 

Let us now consider the more physical limit of low temperature when f3 » 1. Following the same procedure as 
above, we find for the central thermodynamic quantities up to second order in density, first order in kT and after 
restoring dimensions 


^0 


z = e '“tA p 


P = pkT 


S = Nk 


1-f 


1-h 




SkT'^ ^ X^p 15kT\ 


8Eo J 2y/2 


16^0 yj ’ 


1-H 


9kT 


16Eo J 




AV 




9kT \ 


(30) 


Here A = is the standard thermal length. Note that the higher order corrections in the density consists of 

ordinary quantum corrections that are also present in the commutative gas as well as non-commutative corrections 
that exhibit an explicit dependence on the non-commutative energy scale Eq. Also note that the dimensionless 
chemical potential is given by 


p = i log 2 = -1 -I- i log (A^p) -I- log 


^ 3 \ A^p / 15 

8/3 J 2 V 2 V 16/3 


(31) 


The above approximation is only valid when p, < —1, which places constraints on the value of A^p for the validity of 
the approximation. In particular note that one cannot take the zero temperature limit before the zero density limit. 
Finally we consider the case when 0 > p = i log 2 > — 1 and low temperature ,0 >> 1. The way to proceed in this 

case is to perform a Sommerfeld expansion, which we implement here in the following form [^: 


J dxf{x)\og{l + e ^^)=(3 J dxf{x){p-x) + ^f{p) + ... 


(32) 


This yields 


2M3 / TTy/l - fP 
Stt y ey 


1 

12 


^2-\/l - p2 _|_ 2 ) -p 37rp + 6p 



(33) 


To solve for the chemical potential we write pasp = —l-|-po + ^ and expand the g-potential in orders of /3 ^ and 
po up to second order, which assumes a low density. This yields 


p, = —I + ep 





(34) 
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where ep = Egep = is the Fermi energy. This expansion is only valid if ejp > (l — ^). The 

thermodynamics now follows as usual and after restoring dimensions we find 


P = 

kNn'^ 


1 + 7r2(A:T)- 


/ 5 


Vl2e2, ASEoep 


29 \ 

896 E^ J 


2/3eF 


1 - 


£f 

4i?n 


32E§ 


(35) 


Finally we remark that all the results above agree with the results for a commutative perfect Fermi gas [5] when 
0 —)■ 0 and Eg —>■ oo. 


2. The high density limit. 


Now that we have established the low density thermodynamics, we can construct the thermodynamic quantities at 
high density straightforwardly through the use of the low/high density duality. In particular we note from (271 


_^dq{Mj,z-\u!) M3 


— Z- 


dq{ M,P,z,ui) 
dz 


m3 


(36) 


Here Nh is the number of holes or unfilled states. This implies that the thermodynamics at high density is determined 
by the thermodynamics of a gas of holes at low density and in particular the fugacity of the low density gas of holes 
is related to the inverse fugacity of the dense gas of particles, i.e., zt = z~^. The essential difference between the 
high and low density limits resides in the logarithmic term of the fugacity where the inverse fugacity is determined 
by Nh through (36). Since we have already computed the particle fugacity at low density, we can simply read off the 


corresponding inverse or hole fugacity at low density by replacing N by Nh and correspondingly p, p hy ph = ^ and 

Ph = ^, respectively. Also note the following relations between the dimensionless and dimensionful particle and hole 
densities 


p + Ph = l, 

P T Ph — Pmax- (^'^) 


The dimensionful pressure and entropy of the high density gas can now be expressed quite simply in terms of the 
corresponding quantities of the low density hole gas. In particular we have from (271 and the thermodynamic relations 

(pa 


P(M,^, Z,W) = -pmaxfcrlogZ/, +P(M,^, Z/,,£l)), 

S{M,P,z,Cj) = S{M,/3,Zh,uj). (38) 


From these results it is simple to write down the high density thermodynamics dual to the low density cases studied 
above. One simply applies the results of the low density case with the replacements N ^ Nh and p ^ ph- The only 
additional term, which contains the essential physics at high density, appears in the pressure as a logarithm. 

Let us consider the high temperature and high density case where (7=4 log z > 1. As the duality z -H- z“^ translates 


on the level of the chemical potential io fl —p, this implies for the dual low density gas of holes P = j logz < — 1. 
We can then use the results of (29) and (38) to write down the hole fugacity, pressure and entropy explicitly. However, 
the detailed expressions are not relevant for our current discussion and we focus here only on the dominant terms 
that highlight the effects of non-commutativity. Other terms, corresponding to the contribution from the low density 
gas of holes and correction terms that arise from expanding the logarithm to second order in the hole density, are 
represented by dots in the expressions that follow. We find 


P(M,j3,p,Cb) = -pniaxfcriog 



(39) 


Here the leading logarithmic term has been written explicitly in terms of the particle density using (37). We note that 


the pressure diverges logarithmically at the maximum density. It is also insightful to make the volume dependence of 
the pressure explicit at fixed particle number by noting that 

Pmax » 


P{M,j3,V,Uj) = -pma.xkTlog 





(40) 


V 
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Again the divergence in the pressure is evident as the minimum volume is approached. This implies that these 
systems cannot be compressed beyond a minimum size, which corresponds to the minimum number of states required 
to accommodate the number of particles. 


Similarly we find in the case of high density A = ^ log z > 1 and low temperature from (30) and (381 


/5, p, Ca) — PmaxAg Pniax^lf^log (A (pniax p)) T . . . 


(41) 


Again the pressure exhibits a logarithmic divergence as the density approaches the maximum density. It is also 
important to note the large temperature independent degeneracy pressure that contributes even at zero temperature. 
The volume dependence at fixed particle number can also be made explicit as above: 


P(M, P, V, oj) = PmaxA^O - log (A^Pmax) “ Pmax^^log ( 1 - 




min 


+ ... 


(42) 


Once again the pressure diverges at the minimum volume, exhibiting the incompressibility of the system as the 
minimum size is approached. 

Finally we consider the low temperature and high density case when 0 < p = i log 2 < 1, which implies for the 
dual low density gas of holes 0 > p = i logz > —1. Recognising the logarithmic term in the duality relation as 


/3p = —/3p;i, the pressure and entropy can be written down immediately by collecting the results from (34) and (35). 
We only show the logarithmic contribution to the pressure explicitly 


PiMj,V,Q)=p„ 


En-e1 


7r2(fcT)2 

12e^ 


1 - 


2Eo 


+ 


(43) 


Here = (f^)^/^^ is the Fermi energy of holes. We note again the presence of a temperature independent 
degeneracy pressure. Some care is required to investigate the behaviour of the pressure when the density tends to the 
maximum density or, equivalently, when the hole density and thus the hole Fermi energy tends to zero. One cannot 
take this limit before the zero temperature limit as the validity of the approximations made here is then violated. 
However, one can take the zero temperature limit followed by the zero hole density limit, in which case the pressure 
shows no divergence. 

The high/low density duality has important repercussions for the entropy at high density. In particular we note 
from ([M|), (|2^ and (30) that the entropy at high density is extensive in the number of holes, iV/j, rather than the 


number of particles. This implies that the entropy at high density generally exhibits a non-extensive behaviour in 
the particle number, mass and volume of the system. The precise dependence of the entropy on these quantities 
will generally depend on the physical process under consideration. One such process is one in which the system is 
compressed by an external pressure to a finite number, n, of radial units (recall that the radius is quantised as in 
@) larger than the minimum radius M„iin (see (H). Particles are now added to the system leading to an increased 
minimum size, but the size of the system is maintained at the same number of radial units above the increasing 
minimum size by adjusting the external pressure. Under these conditions the number of holes is given by 


m3 

Nh = — -N = V-V^i^ = 


(M, 


-nf 




= nMl^ + 0(M-ij = 


nA„ 


(44) 


3 3 --mm . - ' 

Here is the dimensionful area of the minimal including volume. For this process, this is therefore the quantity 
on which the entropy depends extensively at high density and thus the entropy will exhibit an area law rather than 
a volume dependence at high density. 

Here our interest is in the behaviour at low temperatures. In addition we are interested in the extremal high density 
case where the hole density is low, which corresponds to a chemical potential jl > 1. Then (30) applies to the gas of 
holes and from (38) we can write the entropy to leading order in Amin 


5 = 


5nk 


(45) 


In section we consider a gas confined by gravity, which results in a different non-extensive behaviour for the 
entropy. 


B. Numerical results 

In this section we present exact numerical computations for the thermodynamic quantities described in the pre¬ 
vious section, compare with the various approximations and point out the generic features of these thermodynamic 
quantities, with emphasis on the non-commutative effects. 
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Figures ([^ and ([^ summarise our results. All quantities shown in these graphs are dimensionless. Figures ([^) 
and &) show the dependence of the pressure on the density at low and high temperatures, respectively. Shown on 
these graphs are also the approximations based on equations (29) and (30) and their high density duals. Overall 
there is good agreement between the exact and approximate results in their region of validity. Figures &) and ([^) 
display the same information, but now for the temperature dependence of the pressure. Figure ([^) demonstrates the 
saturation of pressure at high density and low temperature at the temperature independent dimensionless degeneracy 
pressure P = 1. Figure ([^) shows the dependence of the pressure on volume for three different temperature values. 
Here a dimensionless minimum volume of Hmin = 10 was chosen. In all cases the incompressibility at minimum system 
size is apparent. Also note that these curves represent the isotherms. 

Figures (Isk) and dsl: 


show the dependence of the entropy on the density at hxed volume V. The quantity shown 
in these graphs is ^ . Also shown are the low temperature and high temperature entropies at low and high densities 
as computed in equations (29) and (30) and their high density duals. Figures ([^) and ([^) show the same quantities 
at low and high densities as a function of the dimensionless temperature. There is good agreement in the region of 
validity of the approximations. 

It is useful to also study the response of such a gas under isothermal and adiabatic compression. Figure ([^) shows 
the isothermals, i..e., the pressure-volume relation when the gas is compressed at fixed temperature. It exhibits the 
normal behaviour, except for the existence of the minimal volume where the pressure diverges. Figure ([^) shows 
the volume-temperature relation under adiabatic compression. Again the behaviour is standard, except that the 
volume is again bounded from below by the minimal volume. Figure ([^) shows the pressure-volume relation under 
adiabatic compression. This is the most interesting result as one observes a divergence of the pressure at every value 
of entropy. This stems from the fact that entropy decreases with decreasing volume and eventualy vanishes at the 
minimal volume. However, if the entropy is fixed, i.e., no heat is allowed to flow into or out of the gas, the system 
cannot reach volumes with entropy below the initial entropy, which limits the finite volume and manifest itself as a 
diverging pressure. 


C. The ensemble with Ll°^ 7 ^ 0. 


The study of the ensemble with ^ 3 °* 7 ^: 0 is much more involved as the a and A integrals can no longer be performed 
explicitly. However, much progress can be made by realising that the rotating ensemble can be mapped onto an 
unconstrained ensemble with an effective density of states. This effective density of states can then be used to analyse 
the behaviour of these ensembles. We therefore begin by deriving the equivalent unconstrained ensemble and effective 
density of states. The starting point is to rewrite the g-potential (19) in the following way 




dX 


da 


dx 9{l — X^ — x^)0(A^ — a^)0(A) d{x, A) log 


1 + 


(46) 


This suggests that a transformation from the variables x,X,a to new variables x'^X',a' can be found that leaves 
the domain of integration invariant, while transforming away the a dependence in the exponential appearing in the 
logarithm. It is simple to see that the desired transformation must satisfy the following conditions: 


X — (ha = \/l + h^x ', 
-I- A^ = x'^ -I- A'^, 

A2 = A'^ -a'l 


(47) 


Solving these conditions yields 


x{x', A', a') = 




a{x', A', a') = 




\lh (wx'^ — 2a'x' — ha'^) + X'^ (h^ + I) 

X{x',X',a') = ^ -. 

+ I 


( 48 ) 
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FIG. 2: Summary of results for dimensionless pressure. Figure (a) shows the dependence of pressure on the density (solid 
line) at low temperature (5 = 2. Also shown are the low density (dashed line) and its dual high density (dash-dotted line) 
approximations at low temperature, (b) Shows the same result (solid line) at high temperature fi = 0.1. (c) Shows the 

temperature dependence of the pressure (solid line) at low density. Also shown are the high (dashed line) and low temperature, 
low density (dash-dotted line) approximations, (d) Shows the same result at high density, (e) Shows the exactly computed 
pressure at high density (p = 0.99). One notes that the pressure saturates at the degenerate value of P = 1 at low temperature, 
(f) Shows the volume dependence of the pressure at fixed temperatures {P = 0.1 (solid line), P = 1 (dashed line) and P = 2 
(dash-dotted line)) for a minimum system size Knin = 10. 


Conversely 


x\x, A, a) 
a'(x, A, a) 


X — cua 

a + ujx 
•\/w^ -I- 1' 

a/w {Pjx'^ + 2ax — wa^) -I- A^ (w^ -|- 1) 
+ 1 


X'{x, A, a) 


(49) 
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(c) 



(d) 



FIG. 3: Summary of results for dimensionless entropy. The quantity shown is Figure (a) shows this at low temperature 
(/3 = 5). Also shown are the low (dashed line) and high (dash-dotted line) density approximations at low temperature, (b) 
Corresponding result at high temperature (/3 = 0.1). Also shown are the low (dashed line) and high (dash-dotted line) density 
approximations at high temperature. Figures (c) and (d) show the corresponding results at low (p = 0.1) and high densities 
(p = 0.9) as a function of the dimensionless temperature T = i. In figure (c) the high (dashed) and low (dash-dotted) 
temperature approximations at low density and in (d) the high (dashed) and low (dash-dotted) temperature approximations 
at high density are also shown. 


One can check the reality of the solntion for A, X over the whole integration domain and positivity is ensured by the 
choice of positive root. Substituting this transformation into the integral, taking into account the Jacobian, leads to 


where 


/ I py/l—x''^ nX' 

dx' / d\' / da'd'{x', a',C j) log 

-1 Jo J-x’ 


1 -h 


d'{x', X', a', uj) = 


aVi - A'^ - x'^ 




(50) 


(51) 


This is exactly the g-potential of an unconstrained system (w = 0) with inverse temperature /3-\/l + 0^ if we identify 
the effective density of states 


dcs{x',oj)= f dX' f da'd'{x', X', a', ui), 

Jo J-X' 


X' 

-X' 


(52) 


although in this case the X' and a' integrals cannot be explicitly computed (actually the a' integral can be performed, 
but not the A'). It is also useful to note that this can be written as 


dcs{x',uj) = 


VT 


OJ 


( 2,2 rx/i-x''^ ru+ 

— / dX' / dud'{x'j X\u^uj), 

Jo Ju- 


(53) 
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FIG. 4: Response of the system under adiabatic compression. Figure (a) shows the dimensionless volume-temperature rela¬ 
tionship {T = j) and (b) the pressure-volume relationship. 


where and 


d'{x', X', u, Cj) = 


aVi - 

TT (1 — V?) 'JX'^ + 


(54) 


The thermodynamics following from this g-potential can now be studied as before. For reasons of economy we do 
not consider all the possible limits, but only the physically most interesting ones of low temperature (/3 >> 1) and 
slowly rotating gases (|a;| << 1). In this approximation the q-potential reduces to 


q{M,P,z,ui) = 


2M3 

Stt 




(55) 


Before proceeding, we consider the high/low density duality on the level of this approximated g-potential. It is a 
simple matter to show that this q potential can be rewritten as 


q{M,(3,z,Cj) = — ( 1 -b ^ ) log z + q{M,l3,z~'^,oj). 


(56) 


Let us now consider the low density limit with jd = \ logz < — 1. Then the g-potential can be approximated as 


q{M,P,z,uj) = 


2M^z ^ ~,2ro^'2 


Stt 


dx' Vl — x'‘^ ( 1 + 


10(1-x'^) / 


(57) 


Performing the integration and restoring dimensions the following results for the thermodynamics to leading order in 
p and /3 


P = pkT, 


S = Nk 


U = NkT 


325L 




2 - og ( P) - 784£;g,v2moR2 

I 


(tot) 


55L. 


98EoN^moR^ 1 


r T (tot) ^ 

-NEo+ 3 


UNmoR^ ■ 


(58) 
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It is interesting to note that although the pressure does not exhibit any non-commutative corrections due to the linear 
approximation in the density, other thermodynamic quantities do exhibit such corrections for the rotating gas. In 
this case it is also insightful to compute the energy as the rotational contribution can be identified explicitly. In this 
regard also note the shift in the energy by NEq, which simply originates from the shift in energy introduced in (14|. 
Further useful physical information can also be extracted from the expression for the dimensionful u: 


5L 


(tot) 


LJ = 


INmoR^ ’ 


from which we identify the moment of inertia as 


1 = 


TNmoR^ 


(59) 


(60) 


Next we consider the high density dual where M ^ log^ > 1. The g-potential can then be obtained from (56) 
where the low density g-potential is computed as in the case above with M = ^ logz < — 1. One can solve Cj and 2 ; to 

leading order in and p. The expression for 2 is not very informative, but the expression for the dimensional uj 

contains useful physical information as the moment of inertia can be identified from it. For this quantity we find 


3l: 


(tot) 


(jJ = 


^T^PmaxmoR^(l - i \og{\^ph)) 


From this we identify the moment of inertia as 


27r 1 

I = -TTPmax-moR^ {I - ^ \og{X^Ph)) ■ 
O jj 


(61) 


(62) 


Recall that the moment of inertia of a rigid sphere of radius R and matter density p is / = pR^. Comparison with 
(62) suggests that the gas behaves like a body at maximum mass density, yet the deviating factors suggest a form 


of non-rigidity and/or inhomogeneity. Note that the logarithmic term can be neglected at very low temperatures 
and finite hole density. Interestingly, though, the moment of inertia diverges at fixed temperature and vanishing hole 
density. 

For the other thermodynamic quantities we find, after restoring dimensions and keeping only the lowest terms in /3 
and ph, 


P = Eop 

max + kTph - kTp max \og{X^Ph) - 


3kT 


f = iV»(j-log(AV)) + 


SEoTTmoNhR^il - j log(X^ph)) 


2^3 


(tot) 


3 (3 - 2 log {X^Ph}) 


U = -NhkT - NhEo - - - 

2 OTTTOoP; 


16TrEomoPrr, 

3 


llog(A3p,))2 3 

- (tot)^ 


- (tot)"^ 


(63) 


llog(A3p,))2 3 

Note that, as before, one has to be careful with the order of limits here, i.e., one cannot take the T —>■ 0 before ph —t 0. 


V. NON-COMMUTATIVE GAS CONFINED BY GRAVITY 

From the previous sections it is clear that the behaviour of a non-commutative gas at low temperatures and densities 
does not deviate significantly from that of a commutative gas, which makes any observation of non-commutative effects 
extremely difficult. However, at high densities and low temperatures the behaviour of a non-commutative gas deviates 
significantly from its commutative counterpart, which must have observational consequences in systems where these 
conditions can be realised. This may be the case in very dense astro-physical objects, such as white dwarfs and neutron 
stars, although it is unlikely that the densities of these objects are even high enough for non-commutative effects to 
manifest themselves. Despite this, it would still be interesting to study how the non-commutative, non-relativistic 
equation of state alters the behaviour of a gas confined by gravity. Of course, one must realise that such a calculation 
is not realistic in that relativistic effects are important in the study of these systems. Yet, one can expect to gain 
some insight into the qualitative features of these systems through a non-relativistic calculation. 
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Motivated by these considerations, we consider here a non-commutative, non-rotating and non-relativistic gas, 
confined in a large volume, in the presence of gravitational interactions. One expects that the gas will not fill the 
entire volume, but that gravity will confine it to a smaller volume. Our interest will be in the physical properties of 
this confined gas, particularly the density distribution and the mass-radius relation. 

The starting point of this calculation is, as usual, the condition of hydrostatic equilibrium, which reads in the 
Newtonian limit, appropriate for a non-relativistic computation. 


dP{r) 

dr 

dM{r) 

dr 


GM{r)p^{r) 

j,2 

Airr'^Pmir). 


(64) 


Here P{r) is the pressure, M{r) the included mass, Pm{f) the mass density at radius r and G the gravitational 
constant. The mass density is related to the particle density by Pm{f) = Wop(r) where uiq is the mass of the 
constituent particles. It is useful to first rewrite this in dimensionless form. Int rodu cing the dimensionless radius 


r = 5, the dimensionless pressure, P, and particle density, p, introduced in section 


III 


this can be cast in the form 


dP(f) 

dr 

dM{f) 

dr 


GM{r)pir) 

~ f2 

f^p{r). 


(65) 


Here M = ^ and G = 


Gmf,9 


are the dimensionless included mass and gravitational constant. 


mo. ~~ h? 

If the Fermi energy is large compared to A:T, which is usually the case in dense astro-physical objects, we can ignore 
temperature variations and treat the temperature as approximately zero and independent of f. As the density has a 
r dependence, this must be inherited by the dimensionless chemical potential, i.e., p.{f). From (28) it is simple to see 
that 


which implies from (651 


dP{f) djl{f) 
dr dr 


M{f) 


r^ dp,{f) 
G dr 


From this it is easy to derive an equation for fl{r): 

d'^pjf) 2 dp{r) ~ 

df2 f pif 


where p{f) is explicitly given by 


p{f) 


1 r . 

TT 1 _|_ g/3(a;-/i(r-)) 


The boundary conditions are 


m ( 0 ) = Po, 


dp{f) 


dr 


= 0 . 


r =0 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 


The value of po determines the central density and the second condition follows from the requirement of a vanishing 
pressure gradient at the origin and equation (66). 

This equation can only be solved numerically and figure shows the mass-radius relation for two dimensionless 
temperatures. Here the radius has been defined as the point where the dimensionless chemical potential reaches the 
value p = —1 and the density and pressure become exponentially small. Correspondingly the mass is the mass included 
within this radius. In contrast to the commutative case, where the hydrostatic equilibrium condition only admits a 
solution below a threshold mass (the Chandrasekhar mass), there is no threshold mass in the non-commutative case. 
This is due to the incompressible behaviour at maximum density, already alluded to in section [TV] when all available 
single particle states are occupied and the pressure diverges, thereby stabilising the system against collapse. This can 
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R 


FIG. 5: Mass-radius relationship for /? = 10 (dashed line) and P = 50 (solid line). For demonstrative purposes the dimensionless 
gravitational constant was taken as (the unrealistically large) value G — 10“®. All quantities shown are dimensionless and the 
quantity shown on the vertical axis is GM. 

(a) (b) 




FIG. 6: (a) Density profile for po = 5 (solid line), po = 10 (dashed line) and po = 15 (dash-dotted line) for G = 10 (b) 

Entropy densities for the same parameter sets as in (a). All quantities shown are dimensionless. 


only happen when the dimensionless central chemical potential and pressure exceed the critical value of one and the 
dimensionless central density saturates at the maximum value of one. In this case one expects the system to behave 
as an incompressible liquid drop. 

To verify that this is indeed the case we show in figure]^ a) the density distribution for three values of central chemical 
potentials larger than one. From this figure it should be clear that the confined gas behaves as an incompressible 
liquid drop with maximum density in the bulk and a rapidly vanishing density at the edge. Indeed, one notes that 
the edge becomes thinner with increasing radius, indicating that at large mass and radius the incompressible fluid is 
surrounded by a very thin dilute gas of particles. This is very reminiscent of quantum Hall liquids at incompressibility, 
which are known to be closely related to two dimensional non-commutative systems. The situation encountered here 
therefore seems like a natural three dimensional generalisation of quantum Hall fluids. 

Due to the incompressible nature of the bulk, one expects that the excitations of such a fluid will be described by 
the edge excitations. This is certainly the case for quantum Hall fluids. One way of verifying this is to define an 
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FIG. 7: (a) Log-log plot of mass against radius for /3 = 50. The dotted lines both have slopes of three. This demonstrates 
that the dimensionless mass is extensive (M ~ V) in both the low and high density regions, (b) Log-log plot of dimensionless 
entropy S = S/k against dimensionless mass. The left dotted line has a slope of one and the right dotted line of 1/3. This 
demonstrates different scaling behaviours of the entropy in the low {S ~ V) and high (S ~ density regions with a cross 

over between the two. 


entropy density as 


s{f) = i f,(f) - 

V\ d/3 


(71) 


where q{r) is the g—potential of (28) evaluated at /i(r). This is shown in figure [^b). From this it is clear that the 
entropy vanishes in the bulk and is localised in the thin shell on the surface of the droplet. This suggests that this is 
also where the active degrees of freedom of the system must reside, as one would have expected. 

The evolution of the system with increasing mass is further clarified in figure]^ Figure]^ a) shows a log-log plot 
of dimensionless mass against dimensionless radius. The two dashed lines both have a slope of three. This clearly 
shows that in the low and high density regions the dimensionless mass scales with the dimensionless volume, but with 
different scaling factors, i.e. densities, as demonstrated by the different intercepts of the two lines on the vertical axis. 
These two regions correspond to a low density gas and an incompressible liquid drop, as explained above, while there 
is a cross-over region (the s-curve) in-between. This cross-over is associated with a change in the scaling behaviour of 
the entropy, defined as the volume integral of the entropy density (71) over the system volume. This is demonstrated 
in figure 0b). which shows a log-log plot of the dimensionless entropy S = S/k and dimensionless mass. The left 
dotted line has a slope of one and the right one a slope of 1/3. This implies S V, i.e. normal extensive behaviour, 
in the low density region. In the high density regions, however, the scaling is linear in the radius, i.e S ~ and 
non-extensive. This is somewhat counterintuitive in that one would have expected a quadratic dependence or area 
law. In this regard it should, however, be pointed out that the entropy calculation done here is somewhat naive in two 
respects: The first is the computation of the entropy as an integral of a rather naively defined entropy density 
A more appropriate approach, as used in quantum Hall liquids, would be to derive the effective theory describing the 
surface excitations of the incompressible fluid and compute the entropy from there. The second is that the current 
computation completely ignores the quantised nature of the radius, which should become important in the case of 
very thin shells with thickness of the order of a few radial quanta. At this stage it is unclear how these issues should 
be addressed and they will be pursued elsewhere. 


VI. CONCLUSION 

We have studied the thermodynamics of a non-rotating and slowly rotating gas of fermions confined to a spherical 
well in three dimensional fuzzy space, with emphasis on non-commutative effects. These corrections were computed 
at low temperature and density where the non-relativistic approximation should be perfectly valid. These corrections 
may provide useful clues in attempts to observe the effects of non-commutativity at low energy scales and low densities. 

As expected non-commutative effects are manifest at high temperatures and densities. The most prominent of 
these are the existence of a minimal volume at which the gas becomes incompressible and the fact that the entropy 
may exhibit non-extensive behaviour at high densities. These systems also exhibit a remarkable low/high density 
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duality, which underpins many of the features observed at high density. This duality, in turn, is a manifestation of 
an infra-red/ultra violet duality in the single particle spectrum. 

Finally a gas confined by gravity was investigated. The most remarkable feature that emerged in this case is the 
formation of an incompressible liquid drop at high central pressure. This droplet is saturated at maximum density in 
the bulk with a rapid decay to a low density at the edge, giving rise to a very thin, dilute gas of particles surrounding 
the droplet. Entropy calculations suggest that the latter carries the active degrees of freedom of the droplet, which 
gives rise to a non-extensive behaviour of the entropy as a function of the droplet size. 

This study focussed on infinite wells. Finite wells and other confining potentials offer even more intriguing possibil¬ 
ities. In particular, finite wells of sufficient depth have a gap between the finite number of bound states and scattering 
states, which also terminate at high energies [7j. The scattering states can therefore act as a reservoir for the finite 
number of bound states. This interplay may have interesting dynamical consequences that will be studied elsewhere. 
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